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Abstract 

With the help of a phenomenological approach outlined in the text in some detail, we have 
dealt here with the description of the plots on rapidity and pseudorapidity spectra of some 
hadron-secondaries produced in various nucleus-nucleus interactions at high energies. The agree- 
ment between the measured data and the attempted fits are, on the whole, modestly satisfactory 
excepting a very narrow central region in the vicinity of y=rj=0. At last, hints to how the steps 
suggested in the main body of the text to proceed with the description of the measured data 
given in the plots could lead finally to a somewhat systematic methodology have also been made. 
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1 Introduction and Background 

In a chain of our previous works we studied extensively the properties of the rapidity (pseudora- 
pidity) spectra of the various secondary particles in some very high energy collisions. The left- 
overs of the available latest data on rapidity (pseudorapidity) spectra would here be dealt with 
for some specific secondaries produced in some very high energy nuclear collisions. The secondary 
species studied herein are mostly clearly identified particles. Our objective here is quite clear. In 
the light of a Grand Combinational Model (GCM) we would dwell here upon some aspects of the 
behaviour of the rapidity (pseudorapidity) spectra of the identified secondaries produced in high 
energy Au+Au and Pb+Pb interactions. This would provide us an opportunity to check up the 
role of Grand Combinational Model (GCM) in explaining some very recent and interesting data on 
Au+Au and Pb+Pb reactions. 

Amidst the two previous works, this paper is more in line with our second workpQ on under- 
standing the nature of the rapidity spectra for production of some heavy baryons than with the 
first one. In this paper, an empirical energy-dependence of one of the parameters was introduced 
and a systematic approach to the study was built up. This work is just a follow-up of that par- 
ticular methodology. This work is essentially just an exhaustive study on the rapidity-density or 
pseudorapidity-density of the identified hadronic secondaries produced in some high energy nuclear 
collisions. This apart, some clues to developing this procedure as a systematic approach have also 
been highlighted. 

The paper is organized as follows. In the next section (Sec. 2) we give the basic outlook and the 
approach to be taken up for this study. The following section (Sec. 3) provides description of the 
data analyses on Au+Au and Pb+Pb interactions mostly in the graphical plots. The last section 
is reserved for summing up the conclusions with some suggestive remarks to develop the applied 
procedure into a somewhat complete systematic methodology. 

2 The Phenomenological Setting : Premises and the Pathway 

Following FaesslerjJ], Peitzmann[3j, Schmidt and Schukraftjl] and finally Thome et al[5], we [SI [7] 
had formulated in the past a final working expression for rapidity distributions in proton-proton col- 
lisions at ISR (Intersecting Storage Rings) ranges of energy-values by the following three-parameter 
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parametrization, viz, 

a^ = Cl(1 + eXp ^ ) (1) 
where Ci is a normalization constant and yo, A are two parameters. The choice of the above form 
made by Thome et al[5] was intended to describe conveniently the central plateau and the fall-off 
in the fragmentation region by means of the parameters yo and A respectively. Besides, this was 
based on the concept of both limiting fragmentation and the Feynman Scaling hypothesis. For all 
five energies in PP collisions the value of A was obtained to be ~ 0.55 for pions[6] and kaons[7], 
~ 0.35 for protons/antiprotons[7J, and ~ 0.70 for A, S, (ft, £ and f2. And these values of A are 
generally assumed to remain the same in the ISR ranges of energy. Still, for very high energies, 
and for direct fragmentation processes which are quite feasible in very high energy heavy nucleus- 
nucleus collisions, such parameter values do change somewhat prominently, though in most cases 
with marginal high energies, we have treated them as nearly constant. 

Now, the fits for the rapidity (pseudorapidity) spectra for non-pion secondaries produced in 
the PP reactions at various energies are phenomenologically obtained by De and Bhattacharyya[7] 
through the making of suitable choices of C\ and yo- It is observed that for most of the secondaries 
the values of yo do not remain exactly constant and show up some degree of species-dependence . 
However, for A, S, S, 0, and (ft, it gradually increases with energies and the energy-dependence of 
yo is empirically proposed to be expressed by the following relationship |6j : 

y = kin y/sNN + 0.8 (2) 

The nature of energy-dependence of yo is shown in the adjoining figure (Fig.l). Admittedly, as k is 
assumed to vary very slowly with c. m. energy, the parameter yo is not exactly linearly correlated to 
In y/spjN, especially in the relatively low energy region. And this is clearly manifested in Fig.l. This 
variation with energy in k- values is introduced in order to accommodate and describe the symmetry 
in the plots on the rapidity spectra around mid-rapidity. This is just phenomenologically observed 
by us, though we cannot readily provide any physical justification for such perception and/or 
observation. And the energy-dependence of yo is studied here just for gaining insights in their nature 
and for purposes of extrapolation to the various higher energies (in the centre of mass frame, v / *iViv) 
for several nucleon-nucleon, nucleon-nucleus and nucleus-nucleus collisions. The specific energy (in 
the cm. system, ^/snn) for every nucleon-nucleus or nucleus-nucleus collision is first worked out 
by converting the laboratory energy value(s) in the required cm. frame energy value(s). Thereafter 
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the value of yo to be used for computations of inclusive cross-sections of nucleon-nucleon collisions 
at particular energies of interactions is extracted from Eq. (2) for corresponding obtained energies. 
This procedural step is followed for calculating the rapidity (pseudorapidity)-spectra for not only 
the pions produced in nucleon- nucleus and nucleus- nucleus collisions [6J. However, for the studies 
on the rapidity-spectra of the non-pion secondaries produced in the same reactions one does always 
neither have the opportunity to take recourse to such a systematic step, nor could they actually 
resort to this rigorous procedure, due to the lack of necessary and systematic data on them. 

Our next step is to explore the nature of f(y) which is envisaged to be given generally by a 
polynomial form noted below : 

f(y) = a + /3y + 7 y 2 , (3) 

where a, j3 and 7 are the coefficients to be chosen separately for each AB collisions (and also for 
AA collisions when the projectile and the target are same). Besides, some other points are to be 
made here. The suggested choice of form in expression (3) is not altogether fortuitous. In fact, 
we got the clue from one of the previous work by one of the authors (SB) [8] here pertaining to 
the studies on the behavior of the EMC effect related to the lepto-nuclear collisions. In the recent 
past Hwa et al[9] also made use of this sort of relation in a somewhat different context. Now let 
us revert to our original discussion and to the final working formula for ^ in various AB (or AA) 
collisions given by the following relation : 

^Wqx = C 2 {ABY + ^ y2 ^\pp^ Q x = C 3 (AB)^y\l + exp V -^)-\ (4) 
ay ay A 

where C% is the normalization constant and C^=C2{AB) a is another constant as a is also a constant 
for a specific collision at a specific energy. The parameter values for different nucleus-nucleus 
collisions are given in the Tables (Table2 - Table 11). 

However, it is to be noted that the relationship between rapidity and pseudorapidity is given by 
the following standard relation 



dN I m 2 dN 

1 ~ ^2 .2.. A..AJ2. ( 5 ) 



dr\d"p\ V m 2 -.cosh?y dydp, 
with the following properties : 

(a) In the region y»0, %L « ^ 

(b) But, in the region y— >0, there ia a small depression of the t distribution relative to ^ due 
to the above transformation. In experiments at high energies where ^ has a plateau shape, this 



transformation gives a small dip in ^ around i] «0. 

(c) In the cm. frame, the peak of the distribution is located around y~ rj ~0, and the peak 
value of is smaller than the peak value of And this Diminutive Fraction Factor (DFF) is 
given by 

3 Depicting the Results Obtained 

3.1 A Few Pointed Steps 

The procedural steps for arriving at the results could be summed up as follows : 

(i) We assume that the inclusive cross section (I.C.) of any particle in a nucleus-nucleus (AB) 
collision can be obtained from the production of the same in nucleon-nucleon collisions by multi- 
plying the inclusive cross-section (I.C.) by a product of the atomic numbers of each of the colliding 
nuclei raised to a particular function, which is initially unspecified[10j. 

(ii) Secondly, we accept the property of factorization[lJ of that particular function which helps 
us to perform the integral over px in a relatively simpler manner. 

(iii) Thirdly, we assume a particular 3-parameter form for the pp cross section with the param- 
eters C\, yo and A. 

(iv) Finally, we accept the ansatz that the function f(y) can be modeled by a quadratic function 
with the parameters a, j3 and 7. 

3.2 Final Results Delivered 

The results are shown here by the graphical plots with the accompanying tables for the parameter 
values. Here we draw the rapidity-density of pion(-7r), kaon(X), proton-antiproton(A), eft, Q, E, 
A, E for symmetric Pb+Pb and Au+Au collisions and pseudorapidity-density of charged-particle 
(mainly tt + ) for symmetric Au+Au collisions at several energies which have been appropriately 
labeled at the top right corner. In this context some comments are in order. Though the figures 
represents the case for production of pion(7r), kaon(iT), proton-antiproton(A), cfi, Q, E, A, 5, we do 
not anticipate and/or expect any strong charge-dependence of the results. Besides, the solid curves 
in all cases-almost without any exception-demonstrate our GCM-based results. Secondly, the data 
on rapidity(pseudorapidty)-spectra for some high-energy collisions are, at times, available for both 



5 



positive and negative y(ry)-values. This gives rise to a problem in our method. It is evident here in 
this work that we are concerned with only symmetric collisions wherein the colliding nuclei must 
be identical. But in our expression (4) the coefficient f3 multiplies a term which is proportional to 
y and so is not symmetric under y— >(-y). In order to overcome this difficulty we would introduce 
here /3=0 for all the graphical plots (except Fig.12). These plots are represented by Fig.2 to Fig.ll 
for 7T, K, cf>, N, S, S, A, ft in Pb+Pb and Au+Au collision under different conditions. The parameter 
values in this particular case are presented in tables (Table 2 - Table 10). The graphical plots 
shown in fig.2 and Fig. 3 (for /3=0) are for production of 7r~, K~, K + , 4> in Pb+Pb collisions at 
20A GeV, 30A GeV, 40A GeV, 80A GeV respectively. The diagrams shown in Fig. 4 represent 
the production of tt,N,K,T,,A,'E,Q, in Au+Au interaction at S(NN)=7 GeV (for /3=0). And 
the plots depicted for pseudorapidity-spectra in Fig. 5 to Fig.ll are based on the production of 
charged particle (we consider only 7r + ) in Au+Au collision for different centrality bins at 19.6 GeV, 
62.4 GeV, 130 GeV and 200 GeV respectively (for (3=0). The plots shown in Fig.12 are for the 
production of charged particle (mainly ir + ) for four different energies i.e., 19.6 GeV, 62.4 GeV, 
130 GeV, 200 GeV respectively and the parameter values are shown in Table 11. Here we would 
mention that the data are for 19.6 GeV, 130 GeV, 200 Gev for PHOBOS Collaboration and that 
of 62.4 GeV for STAR Collaboration as shown on the top right corner of the figure. Finally, the 
diagram in Fig. 13 represents the variation of f3 and 7 with the energy values and we draw a mean 
curve in this particular diagram. 

4 Concluding Remarks and Some Comments 

On an overall basis, our model-based results are in fair agreement with the most of the data- 
sets, excepting y«0 or r/ss0 region, wherein the data shows flat-plateau structures in almost all 
the diagrams exhibiting data on both positive and negative rapidities or pseudorapidities. The 
degree of disagreement in the vicinity of rjxiO region is evidently much stronger for the plots on 
pseudorapidity-density versus pseudorapidity plots. These discrepancies might probably be ascribed 
to our simplistic assumption of y=rj. Had we been able to compute the diminutive fraction factor 
(DFF) as given by expression (6), we would have been capable of giving the pseudorapidity- figures 
much better looks. And this computation is not possible because of the fact that the rapidity- 
data-sets do not generally offer even the slightest hints on the py-ranges of the secondaries under 
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observations and/or measurements. 

The last figure of this paper carries some special physical significance which we now explain 
below. This is, by essence, undoubtedly a purely phenomenological model with no or very little 
predictive capacity. The energy-dependences studied in Fig. 13 for some of the involved parameters, 
j3 and 7, could provide us some insights into what could be the possible values of f3 and 7 at 
some higher /lower /intermediate values of the cm. energies of the interactions for any specific 
secondary. This could help, we believe, to reduce the elements/components of phenomenology and 
introduce some degree, however low, of predictivity of values of f3 and 7 by necessary intrapolation 
or extrapolation, as the case may be, for any specific secondary produced in the same nuclear 
interactions. Thus, if sufficient and reliable data at, at least, six to seven cm. energies at reasonable 
intervals are available allowing the scopes for studying the nature of cm. energy-dependence 
of these parameters, the present procedure could be nurtured to a better and more competent 
methodical approach. 
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Table 1: Variation of yo with Energy. [Reference Fig. No.l] 



Energy (^>s NN ) (UeV ) 


Constant(k) 




&.Z{20AGeV) 


2.76 


5.894 


7 


2.65 


5.951 


7.6(30AGeV) 


2.54 


6.006 


8.7(40AGeV) 


2.40 


6.085 


12.3(80AGeV) 


2.16 


6.276 


19.6 


1.92 


6.517 


62.4 


1.54 


7.153 


130 


1.39 


7.556 


200 


1.32 


7.794 



Table 2: Values of different parameters for production of identified hadrons in central Pb+Pb 
collisions at E^ eam = 20 AGeV (for /3=0) for both +ve and -ve rapidities. [Reference Fig. No. 2(a) 
& 3(a)] 



Production 


c 3 


7 


xL 

ndf 




84.941 ±0.2108 


-0.044 ±0.0003 


21.509/22 


K+ 


20.574 ±0.1556 


-0.057 ±0.0008 


4.119/05 


K~ 


05.507 ±0.0447 


-0.091 ±0.0015 


7.753/14 




01.325 ±0.0406 


-0.139 ± 0.0021 


4.423/07 



Table 3: Values of different parameters for production of identified hadrons in central Pb±Pb 
collisions at Eb ea m = 30 AGeV (for /3=0) for both ±ve and -ve rapidities. [Reference Fig. No. 2(b) 
& 3(b)] 



Production 


c 3 


7 


xL 

ndf 




97.912 ±0.2754 


-0.037 ±0.0003 


10.898/13 


K+ 


22.849 ±0.1128 


-0.049 ±0.0005 


5.824/08 


K~ 


07.547 ±0.0705 


-0.075 ±0.0011 


4.316/10 




01.384 ±0.0158 


-0.095 ±0.0008 


4.264/06 



Table 4: Values of different parameters for production of identified hadrons in central Pb±Pb 
collisions at Ei, eam = 40 AGeV (for /3=0) for both ±ve and -ve rapidities. [Reference Fig. No. 2(c) 
& 3(c)] 



Production 


c 3 


7 


xL 

ndf 




111. 998± 0.3811 


-0.035 ±0.0003 


8.926/10 


K+ 


023.363 ±0.1275 


-0.039 ±0.0005 


4.621/06 


K~ 


010.605 ±0.0755 


-0.073 ±0.0003 


3.394/09 




001. 185 ±0.0147 


-0.063 ±0.0006 


4.206/09 
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Table 5: Values of different parameters for production of identified hadrons in central Pb+Pb 
collisions at E^ eam = 80 AGeV (for (3=0) for both +ve and -ve rapidities. [Reference Fig. No. 2(d) 
&3(d)] 



Production 


c 3 


7 


xL 

ndf 




147.686 ±0.4383 


-0.027 ±0.0002 


8.659/11 


K+ 


026.426 ±0.1038 


-0.030 ±0.0004 


5.455/08 


K~ 


012.790 ± 0.0763 


-0.041 ±0.0007 


9.002/12 





001.762 ±0.0044 


-0.044 ±0.0002 


0.543/06 



Table 6: Values of different parameters for production of identified hadrons in central Au+Au 
collisions at y / sj\F/v = 7 GeV (for (3=0) for both ±ve and -ve rapidities. [Reference Fig. No. 4] 



Production 


c 3 


7 


JSL 

ndf 


ir 


299.742 ±0.4296 


-0.026 ±0.0001 


8.014/07 


N 


136.702 ± 0.2569 


-0.048 ±0.0005 


5.460/04 


K 


011. 734 ±0.0248 


-0.030 ±0.0001 


0.621/03 


S,A 


015.886 ± 0.0055 


-0.056 ±0.0004 


4.750/05 


n 


000. 017 ± 0.0004 


-0.079 ±0.0048 


8.754/10 




000.628 ± 0.0037 


-0.057 ±0.0010 


9.502/14 



Table 7: Values of different parameters for production of charged-particle (vr + ) in Au±Au collisions 
at y/sNN = 19.6 GeV (for (3=0) for both +ve and -ve pseudo-rapidities. [Reference Fig. No. 5] 



Centrality 


c 3 


7 


xL 

7i df 


03% - 06% 


363.098 ± 1.764 


-0.0157 ±0.00009 


17.359/12 


06% - 10% 


328.094± 1.195 


-0.0148 ±0.00007 


12.873/19 


10% - 15% 


266.638 ±0.967 


-0.0137 ±0.00007 


21.209/18 


15% - 20% 


207.392 ± 1.130 


-0.0125 ±0.00008 


9.160/10 


20% - 25% 


180. 106 ±0.656 


-0.0123 ±0.00007 


15.549/16 


25% - 30% 


150.634 ± 0.456 


-0.0117 ± 0.00006 


14.747/12 


30% - 35% 


115.939 ± 0.436 


-0.0110 ± 0.00007 


16.130/15 


35% - 40% 


087.982 ±0.405 


-0.0098 ±0.00008 


11.447/13 


40% - 45% 


071.793 ±0.249 


-0.0093 ±0.00008 


16.840/10 



Table 8: Values of different parameters for production of charged-particle (vr + ) in Au±Au collisions 
at y / S7VA? : = 62.4 GeV (for (3=0) for both ±ve and -ve pseudo-rapidities. [Reference Fig. No. 6 k, 7] 



Centrality 


c 3 


7 


xL 

ndf 


00% - 03% 


858. 977 ±2.349 


-0.0131 ±0.00002 


3.527/04 


03% - 06% 


805.167± 1.252 


-0.0135 ±0.00008 


9.541/12 


06% - 10% 


690.250 ±4.553 


-0.0126 ±0.00006 


15.843/13 


10% - 15% 


505.765 ±2.056 


-0.0111 ±0.00004 


15.882/15 


15% - 20% 


397.679 ± 1.516 


-0.0101 ±0.00004 


11.527/12 


20% - 25% 


288.758 ± 1.128 


-0.0086 ±0.00005 


9.042/10 


25% - 30% 


223.869 ± 1.019 


-0.0081 ± 0.00004 


8.544/08 


30% - 35% 


179.875 ±0.496 


-0.0076 ±0.00004 


10.873/11 


35% - 40% 


126.863 ±0.661 


-0.0066 ±0.00005 


10.169/10 


40% - 45% 


112. 144 ±0.404 


-0.0068 ±0.00005 


6.154/05 


45% - 50% 


091.802 ±0.290 


-0.0077 ±0.00004 


11.346/09 
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Table 9: Values of different parameters for production of charged-particle (vr + ) in Au+Au collisions 
at y / S7v r /v = 130 GeV (for j3=0) for both +ve and -ve pseudo-rapidities. [Reference Fig. No. 8 & 9] 



Centrality 


c 3 


7 


xL 

ndf 


00% - 03% 


1017.26 ±4.767 


-0.0106 ±0.00006 


8.489/08 


03% - 06% 


859.123 ±4.007 


-0.0096 ±0.00005 


11.504/11 


06% - 10% 


771.332 ±4.214 


-0.0094 ±0.00005 


11.448/11 


10% - 15% 


553.134± 1.895 


-0.0083 ±0.00002 


11.348/07 


15% - 20% 


490.218 ±2.056 


-0.0084 ±0.00004 


13.212/15 


20% - 25% 


384.377 ± 1.812 


-0.0079 ±0.00004 


9.556/16 


25% - 30% 


342.650 ± 1.736 


-0.0081 ± 0.00004 


14.193/13 


30% - 35% 


264.987 ± 1.080 


-0.0077 ±0.00004 


8.711/13 


35% - 40% 


183.701 ±0.757 


-0.0061 ±0.00005 


3.271/06 


40% - 45% 


165.105 ± 0.243 


-0.0068 ± 0.00001 


11.067/08 


45% - 50% 


115.942 ±0.417 


-0.0065 ±0.00003 


8.366/13 



Table 10: Values of different parameters for production of charged-particle (tt + ) in Au±Au collisions 
at y / stvjv = 200 GeV (for /3=0) for both ±ve and -ve pseudo-rapidities. [Reference Fig. No. 10 & 11] 



Centrality 


c 3 


7 


xL 

ndf 


00% - 03% 


1305. 84 ±6.334 


-0.0095 ±0.00004 


8.374/07 


03% - 06% 


1217.02 ±3.923 


-0.0092 ±0.00002 


2.677/07 


06% - 10% 


1085.85 ±5.787 


-0.0089 ±0.00003 


8.494/08 


10% - 15% 


786. 155 ±5.756 


-0.0083 ±0.00004 


10.695/07 


15% - 20% 


671. 544 ±0.816 


-0.0080 ±0.00001 


13.220/08 


20% - 25% 


494.043 ±4.083 


-0.0076 ±0.00005 


8.479/08 


25% - 30% 


371.719 ±2.104 


-0.0069 ±0.00004 


7.826/06 


30% - 35% 


298. 854 ±2.012 


-0.0069 ±0.00005 


10.232/07 


35% - 40% 


298.430 ± 1.067 


-0.0077 ±0.00001 


19.517/16 


40% - 45% 


142. 160 ±0.609 


-0.0060 ±0.00003 


11.112/10 


45% - 50% 


203.956 ±0.874 


-0.0063 ±0.00003 


13.112/08 



Table 11: Values of different parameters for production of charged-particle (tt + ) in Au±Au collisions 
at four different energies for +ve pseudo-rapiditie only. [Reference Fig. No. 12 



Energy(GeV) 


c 3 




7 


xi 

ndf 


19.6 


361.238 ±3.900 


0.015 ±0.00032 


-0.020 ±0.00009 


8.596/25 


62.4 


778. 667 ±8.000 


0.018 ±0.00028 


-0.018 ±0.00008 


3.847/04 


130 


511. 297± 7.890 


0.027 ±0.00054 


-0.012 ±0.00015 


14.445/25 


200 


557.934 ± 7.087 


0.028 ± 0.00044 


-0.011 ±0.00011 


9.578/26 
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Figure 2: Rapidity distributions of identified hadrons in central Pb+Pb collisions at Eb eam = 20, 
30, 40, 80 AGeV for [3=0. The symbols are the experimental data and the data points are taken 
from [11] and the parameter values are taken from Table 2-Table 5. The solid curve provide the 
GCM-based results. 
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(c) (d) 

Figure 3: Rapidity spectra for (f> in central Pb+Pb collisions at Eb ea m = 20, 30, 40, 80 AGeV for 
/3=0. The different experimental points are taken from [11] and the parameter values are taken 
from Table 2- Table 5. The solid curve provide the GCM-based results. 
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Figure 4: Rapidity distributions of identified hadrons in central Au+Au collisions at y / s/v r /v=7 GeV 
for (3=0. The different experimental points are taken from [12] and the parameter values are taken 
from Table 6. The solid curve provide the GCM-based results. 
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Figure 5: Pseudo-rapidity spectra for tt + for nine centrality bins representing 45% of the total 
cross-section for Au+Au collisions at - v /sjv r /v=19.6 GeV for /3=0. The different experimental points 
are taken from |13j and the parameter values are taken from Table 7. The solid curve provide the 
GCM-based results. 
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Figure 6: Plot of vs. r/ for ir + for seven centrality bins representing 45% of the total cross- 
section for Au+Au collisions at ^/snn=62A GeV for /3=0. The different experimental points are 
taken from [13] and the parameter values are taken from Table 8. The solid curve provide the 
GCM-based results. 
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Figure 7: Plot of -jj— vs. rj for 7r + for four centrality bins representing 45% of the total cross-section 
for Au+Au collisions at y / s/\T/v=62.4 GeV for /3=0. The different experimental points are taken 
from [13J and the parameter values are taken from Table 8. The solid curve provide the GCM-based 
results. 
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Figure 8: Pseudo-rapidity spectra for 7r + for five centrality bins representing 45% of the total 
cross-section for Au+Au collisions at y / s/v r /v=130 GeV for /3=0. The different experimental points 
are taken from [13] and the parameter values are taken from Table 9. The solid curve provide the 
GCM-based results. 



17 



Au+Au(7t + ) 130GeV 
25%-30%*3 h-^h 
20%-25% l— •— l 
GCM — 




Au+Au(7t + ) 130GeV 
35%-40%*3 i— ■— i 
30%-35% I— •— i 
GCM - 




10 



(b) 



Au+Aufa""") 130GeV 
45%-50%*3 h- ■— i 
40%-45% h^-h 
GCM 



10 



Figure 9: Pseudo-rapidity spectra for tt + for six centrality bins representing 45% of the total cross- 
section for Au+Au collisions at - v /saT/v=130 GeV for (3=0. The different experimental points are 
taken from [13] and the parameter values are taken from Table 9. The solid curve provide the 
GCM-based results. 
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Figure 10: Plot of vs. rj for ir + for six centrality bins representing 45% of the total cross-section 
for Au+Au collisions at y / J/v r iv=200 GeV for f3=0. The different experimental points are taken from 
|13j and the parameter values are taken from Table 10. The solid curve provide the GCM-based 
results. 
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Figure 11: Plot of vs. rj for 7r + for five centrality bins representing 45% of the total cross-section 
for Au+Au collisions at y / J/v r /v=200 GeV for f3=0. The different experimental points are taken from 
|13j and the parameter values are taken from Table 10. The solid curve provide the GCM-based 
results. 
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Figure 12: Pseudorapidity distributions of charged particles(for n + ) for various cm. energies 
in Au+Au central collisions. Pseudorapidity distributions for 0-6% central Au+Au collisions at 
= 200, 130, and 19.6 GeV are from the PHOBOS experiment and 62.4 GeV is from STAR 
experiment (/? t^O). The different experimental points are taken from [TJ] and the parameter values 
are taken from Table 11 . The solid curve provide the GCM-based results. 
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Figure 13: Variation of j3 and 7 with increasing energy. All values are taken form Table 11. 
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